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Abstract 

A modified gravity (MOG) possesses a light, neutral vector particle called 
a "phion" associated with a vector field <p^, which forms a cold fluid of Bose- 
Einstein condensates before recombination with zero pressure and zero shear 
viscosity. The energy density associated with this Bose-Einstein condensate 
fluid dominates the energy density before recombination and produces a den- 
sity parameter, Q,^ ~ 0.3, that together with the fractional baryon density 
Qb ~ 0.04, and a cosmological constant parameter Q\ ~ 0.7 yields an ap- 
proximate fit to the data for the acoustical oscillations in the CMB power 
spectrum. The quantum phion condensate fluid is abundant well before re- 
combination and can clump and form the primordial structure for galaxies. At 
late times in the expanding universe, in local bound systems such as galaxies 
ordinary baryonic matter dominates the matter density. For galactic systems 
in the present epoch, the modified Newtonian acceleration law determined by 
MOG describes well galaxy rotation curve data and X-ray cluster mass profile 
data. 



e-mail: john.moffat@utoronto.ca 



1 Introduction 

A relativistic modified gravity (MOG) theory |2j has been proposed to explain 
the rotational velocity curves of galaxies and the X-ray data for clusters of galaxies 
with a modified Newtonian acceleration law, without non-baryonic dark matter. A 
fitting routine for galaxy rotation curves has been used to fit a large number of galaxy 
rotational velocity curve data, including low surface brightness (LSB), high surface 
brightness (HSB), dwarf galaxies and elliptical galaxies with both photometric data 
and a two-parameter core model without non-baryonic dark matter [21 El- The fits to 
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the data are remarkably good and for the photometric data only the one parameter, 
the mass-to- light ratio (M/L), is used for the fitting, once two parameters Mo and tq 
are universally fixed for galaxies and dwarf galaxies. The fitting results for galaxies 
are consistent with the data for the Tully-Fisher relation A large sample of 
X-ray mass profile cluster data has also been fitted jSJ. The gravity theory also 
fits the anomalous acceleration data observed during the Pioneer 10-11 spacecraft 
missions ;GI , and is consistent with gravity observations in the solar system and the 
binary pulsar PSR 1913+16. 

The MOG requires that Newton's constant G, the coupling constant u that 
measures the strength of the coupling of a skew field to matter and the mass fi 
of the skew field, vary with distance and time, so that agreement with the solar 
system and the binary pulsar PSR 1913+16 data can be achieved, as well as fits 
to galaxy rotation curve data and galaxy cluster data. The variation of G leads 
to a consistent description of relativistic lensing effects for galaxies without non- 
baryonic dark matter. In ref . U\ , the variation of these constants was based on a 
renormalization group (RG) flow description of quantum gravity theory formulated 
in terms of an effective classical action. Large infrared renormalization effects can 
cause the effective G, u, \x and the cosmological constant A to run with momentum k 
and a cutoff procedure leads to a space and time varying G, u and /i, where fi = l/r 
and r is the effective range of the skew symmetric field. In the MOG theory [2], the 
action contains a contribution that leads to effective field equations that describe 
the variations of G, uj and fi. In principle, we can solve for the complete set of 
field equations and determine the dynamical behavior of all the fields. However, 
in practice we make approximations allowing us to obtain partial solutions to the 
equations, yielding predictions for the various physical systems considered. 

A modified gravity theory must be able to explain the following cosmological 
data and phenomena: 

1. The cosmic microwave background (CMB) data, including the power spectrum 
for both large and small angular scales on the sky; 

2. The formation of galaxies before and after recombination; 

3. Gravitational lensing; 

4. N-body simulations of large-scale galaxy surveys; 

5. The accelerating expansion of the universe. 

The standard explanation for these phenomena is the cold dark matter (CDM) 
scenario and some form of uniform dark energy that begins to dominate in the 
present epoch of the universe. The problem with this scenario is that 96% of the 
universe's matter and energy budget is invisible. It remains undetected in ground- 
based laboratory experiments, and has yet to be given a satisfactory physical in- 
terpretation [8J. An attempt to explain the CMB power spectrum at the surface 
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of last scattering about 370,000 years after the big bang in terms of the visible 
baryon number density [S], leads to a satisfactory fit to the WMAP and Archeops 
data EI] for the first two peaks in the power spectrum at small angles with 

the baryon density fraction flj ~ 0.04 and the cosmological constant density fraction 

~ 0.95. This scenario requires a neutrino mass contribution with m u ~ 1 — 2 
eV. The baryon density parameter Qb ~ 0.04 is in agreement with recent estimates 
for this quantity at big bang nucleosynthesis. The problem with this pure baryon- 
neutrino scenario is that a third putative peak in the angular power spectrum is 
erased by baryon drag. A problem with neutrinos describing dark matter is that 
they become non-relativistic at a late-time epoch in the expanding universe, and 
their fluctuations are erased in the power spectrum. 

The most recent balloon-born Boomerang data (13] indicate the existence of a 
third peak in the power spectrum, which is erased by baryon drag without cold dark 
matter ^U|. If future data confirms that a significant third peak in the angular power 
spectrum exists, then this strongly supports that some form of dark, non-baryonic 
matter should exist at the surface of last scattering. All of this suggests that some 
kind of non-relativistic dark matter should be accounted for in cosmology. On the 
other hand, the MOG fits to the galaxy rotation curve data and the X-ray galaxy 
cluster data [3J E] without non-baryonic dark matter are remarkably good and avoid 
problematic issues related to CDM fits to the galaxy data. 

It is also necessary to provide a satisfactory explanation for the early universe 
growth of structure and the eventual emergence of galaxies and clusters of galax- 
ies. Large computer N-body simulations show that some kind of non-relativistic 
dark matter has to constitute the proto-galaxies about 200 million years after the 
big bang [HJ US] . The data suggest that early universe density fluctuations were 
a Gaussian random field described by a cold dark matter (CDM). The standard 
assumption is that the CDM is made of elementary particles that interact only 
gravitationally. Baryons and neutrinos alone fail to produce a realistic description 
of large-scale galaxies and clusters surveys. 

An important extra-degree of freedom in MOG is a light, electrically uncharged 
vector particle called a "phion". In the early universe at a temperature T < T c , 
where T c is a critical temperature, the phions become a Bose-Einstein condensate 
(BEC) fluid. The phion condensates couple weakly with gravitational strength to or- 
dinary baryonic matter. This cold fluid of phion condensates dominates the density 
of matter at cosmological scales and, because of its clumping due to gravitational 
collapse, allows the formation of structure and galaxies at sub-horizon scales well 
before recombination. We do not postulate the existence of cold dark matter in the 
form of heavy, new stable particles such as supersymmetric WIMPS. 

In the following, we shall calculate the power spectrum for values I 3> 1 obtain- 
ing the positions and heights of the acoustical peaks from cosmological parameters. 
We shall use the analytical formula derived by Mukhanov JHj and earlier work by 
Hu and Sugijama [T7| . Frampton, Ng and Rohm [TH], Weinberg [TH] and Dodel- 
son 20 . The approximate analytical formula yields results for the small angular 
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scale acoustical oscillations that agree with reasonable accuracy with the WMAP, 
Archeops and Boomerang data [HI H21 CHI- Since the phion condensates do not 
couple to radiation and only couple with weak gravitational strength to baryons, 
their density perturbations are not subject to dissipation like photon and neutrino 
density perturbations before recombination. 

For large-scale cosmology, the phion condensate density on the average dominates 
the matter density, while for locally bound galactic systems and clusters of galaxies 
visible baryon density dominates the matter density, and the MOG field equations 
and the quantum BEC equations lead to an effective modification of the Newtonian 
acceleration law for weak fields that yields an excellent fit to the galaxy rotation 
curve data and the X-ray cluster mass profile data IE] ■ 



2 Classical Field Equations and Test Particle Mo- 
tion 

The relativistic gravitational field equations in MOG are given by (in units with 
c=l)0: 

G^ v - S>A + = 8nGT^ v , (1) 

V V B^ + + -V u uB^ = --J", (2) 

o<p^ UJ UJ 

where G^ u = R^ u — ^g^ u R, A is the cosmological constant and V M denotes covariant 
differentiation with respect to the metric g^. We have 

Qtu> = G(v a v a e 9llu - v M v„e), (3) 

where G(x) is the scalar field spacetime dependent gravitational "constant" and 
Q(x) = 1/G(x). Moreover, 

By, V = dyrffy, ~ d V (j)^ (4) 

and denotes a matter current. The potential V(<p) has the form 

V^) = -^ 2 ^ + W^), (5) 

where fi(x) is the effective spacetime dependent mass of the vector field M and W(4>) 
is a 0^ self-interaction potential. The coupling constant u(x) is the scalar spacetime 
dependent coupling of the skew field B^ v to matter. 
The total energy-momentum tensor is given by 

Tfi V = Tmihv + + (6) 

where T M ^ V1 T^^ V and Tg^v denote the energy- momentum contributions from the 
matter fields, the vector field ^(x), and the scalar fields G(x),uj(x) and fi(x), re- 
spectively. 
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From the Bianchi identities 

VuG^ = 0, (7) 
and from the field equations ((TJ), we obtain the conservation law: 

V U T^ + ^VuGT^ - ^VuQ^ = 0. (8) 

Cx 57TG 

The effective gravitational "constant" G(x) satisfies the field equations: 

^ + am +N .^(r + ^), (9) 

where 

N = -3e(-V a GV a G + V(G) \ + G^V a uV a uj - V(u) 

«(iv„,vw M ) + ^ 

and T = g^T^. The scalar field u(x) obeys the field equations 



+— ( n v Q /iV a /i - v(n) ) + — v Q v Q e, (io) 



dV(u) 



V a V a u + ^— ^ + F = 0, (11) 
au; 

where 

f = -eVaGv^ + G^s^s^ + y^V (12) 

The field /i(x) satisfies the equations 

V Q V>+^^ + P = 0, (13) 



where 



eV a GV aj j + -V>V Q /i + w// 2 G 



(14) 



and the last term arises from the \x dependence of V(4>). 

The equation of motion for a test particle of mass m and charge A when u = 
^ren = constant, where u ren is the renormalized value of u, is given by 

where r is the proper time along the test particle trajectory and 

m ar 

For massless photons the test charge A is zero and photons move along null geodesies: 

^ + ^,-- = 0. (17) 
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3 Bose-Einstein Condensate of Phions 



An important prediction of quantum statistical mechanics is that of a phase tran- 
sition in an ideal gas of identical bosons when the thermal de Broglie wavelength 
exceeds the mean spacing between bosons j2lj. The bosons in the lowest energy 
state are stimulated by the presence of other bosons to occupy the same state, re- 
sulting in a macroscopic occupation of a single quantum state that constitutes a 
macroscopic quantum-mechanical system [221 122] • 

Well before recombination, the boson-phion particle can collapse into the lowest 
energy quantum ground state when the temperature of the phion fluid is below a 
critical temperature, T < T c , and form a degenerate Bose-Einstein fluid that does 
not couple to photons and is not dissipated by photon pressure [2H 123 EE1 123 EE] • 
The fluid has zero viscosity and pressure and forms the dominant energy density of 
matter in the universe at large cosmological scales greater than the scales of galaxies 
and clusters of galaxies. The phion Bose-Einstein condensates (BEC) are the dark 
matter contribution to the cosmological density, which dominates the baryon density, 
fl<£ > where denotes the fractional BEC density, at cosmological scales and 
couple weakly with gravitational strength to baryon matter. 

The phions undergo a second-order phase transition through a spontaneous sym- 
metry breaking mechanism below a critical temperature T < T c . This mechanism 
can generate the effective phion mass and the BEC exhibits a broken symmetry 
and a phase coherence. The condensate wave function is non-vanishing and the 
ground state (vacuum) depends coherently on the phases of the phions, such that 
they are spatially correlated throughout the system. For a quantum state of light 
phions, it is possible to form a non-relativistic quantum BEC of phions in a second 
order phase transition. 

The gravitational collapse (clumping) of the phion condensate fluid to form grav- 
itational potential wells before recombination, allows the early formation of struc- 
tures at small horizon scales of the order of proto-galaxy sizes. These structures 
eventually lead to the formation of clusters of galaxies in a hierarchical scheme of 
structure growth. For this to happen, the density of non-relativistic BEC phions 

dominates the clumping fluid with a sub-horizon size of order of the size of a 
proto-galaxy and with a suitable value of the coherence length R c . 

A Jean's scale Aj separates gravitationally stable and gravitationally non-stable 
modes. The time scale for gravitational collapse is given by r grav ~ 1/ ^/G ren p0. 
On the other hand, the time scale for fluid pressure to dynamically respond is 
Tpress ~ A/c s , where A is the size of the condensate in the fluid and c s is the sound 
velocity. In units with ft and c not equal to unity, for quantum non-relativistic 
BEC matter, the velocity v c in the condensate fluid is determined by the de Broglie 
wavelength Xb = h/m^Vc. The fraction of phions in the coherent BEC ground 
state is about 100% before decoupling at a redshift za ~ 1100, and after decoupling 
baryonic structures can begin to form and the quantum Jeans scale Aj quant can 
become the size of proto-galaxies. 
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Three characteristic lengths define the system: the Compton wavelength, A c = 
h/m^c, the average distance between particles, Xd = where denotes the 

number density of phions, and the thermal wavelength A^ = h/p. The thermal de 
Broglie wavelength is 

where ks is Boltzmann's constant. The non-relativistic limit corresponds to Ay 3> 
A c , while the relativistic limit holds for At <C A c . 

For an ideal non-relativistic Bose gas, the number of phions is given by ED] : 

^^'l ^-^M-. ' (19) 

where e^, denotes the phion energy spectrum and z is the fugacity, related to the 
chemical potential /i c h em by z = exp(/x chem /fcs7 1 ). The fugacity z is solved in terms 
of the equation 

27r,„ ,,/o Z" 00 , 1 2 



^ = ^(2m/ 2 de — i +i/r^-' ( 2 °) 

/r jo z 1 exp(e0//c_B-Zj — 1 V 1 — z 



where V is a characteristic volume. We can solve this equation to obtain 

9a/2(z) t 1 z 

where 



n * = ^r + v—> (21) 



We can rewrite (J2TJ) as 



9n{z) = W)l dx z-^M*)-i (22) 

A T ^ = A T ^-^ 3/2 W, (23) 

which implies that {uq^/V > when the temperature T and are such that 

Ajn* > C(3/2) = 2.612..., (24) 

where ^ 3 / 2 (l) = C(3/2) and ((x) is the Riemann-zeta function. A finite number of 

phion bosons occupy the ground state level with = 0. This defines the critical 
temperature for phion condensation: 

2irh 2 n/ 3 

Tc ~ fc B m C 2/3 (3/2)' (25) 



and the critical phion number density 



= (26) 

An 



Thus, for T < T c a large number of phions occupy the ground state = 0. 

We are required in a relativistic cosmology to take into account the expansion 
of the universe. The solution for the case of an ideal boson gas in a Friedmann- 
Robertson- Walker (FRW) universe is complicated, so authors |JTJ 1221 HUH OH EH] 
have considered a boson gas in a closed Einstein universe. In the spatial geometry 
of a 3-sphere, the metric takes the form: 

ds 2 = c 2 dt 2 - a 2 E [dx 2 + sin 2 X (d6 2 + sin 2 Od<f>% (27) 

where x an d run from to In. The radius of the Einstein universe a E is constant. 
The solution for the spin 1 phion discrete energy spectrum is given by (h — c — 1): 



-!(n 2 + m 2 4) 1/2 ( 28) 



a E 

with degeneracy d n = 2(n 2 — 1) where n = 2, 3, .... 

Let us take the mean phion density at large scale galactic distances to be of the 
order of the critical density, p^ ~ 9.52 x 10~ 30 gcm~ 3 . From (|25jl . we get in the 
condensation region 

nt < f;,^^ 3 ' 5 . (29) 



.2tt^C 2/3 (3/2)/ 

We obtain for the phion condensate density that is comparable to the cosmological 
critical density for a critical temperature at recombination T ~ 10 5 K: 

< 6.3 x 10 7 cm" 3 . (30) 

We get from = p^/n^ for the effective phion condensate mass at the epoch of 
recombination: 

m ~ 8.5 x eV/c 2 . (31) 



4 Bose-Einstein Gas and the Ground State 

A Bose-Einstein state was experimentally verified by the Joint Institute for Labo- 
ratory Astrophysics (JILA) in 1995. Anderson et al. [36j produced a condensate of 
spin-polarized 87 RB atoms in a confining magnetic trap. Densities and temperatures 
of order 10 11 cm -3 and tens of micro-Kelvin, respectively, were achieved. 

We shall assume that the phion boson inter-particle interactions are weak, so 
that at low temperatures in cosmology the de Broglie wavelengths of the phions 
are large compared to the inter-phion separation. The phion-phion interactions 
are dominated by elastic S-wave scattering and we can consider for weak fields 
of gravitational strength only two-body collisions. The scattering length can be 
positive, a > 0, for a repulsive interaction and negative, a < 0, for an attractive 
interaction. The phion self-interaction potential is given by 

*7(x-x') =0<J(x-x'). (32) 
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Here, the strength of the interaction is 



AirK a 

m 



(33) 



The Hamiltonian for weakly interacting phions in an external gravitational potential 

n 2 



V^ rav is of the form 



H 



d 3 xip^(x) 



2m, 



grav V 



V>(x) 



+ - 



d 3 x / dVV> t (x)V> t (x')£/(x-x / )V>(x')# 



(34) 



where ip(x) and ip '(x) are the phion field annihilation and creation operators, re- 
spectively. 

A necessary condition for the applicability of a calculation of the BEC gas is 
that a 3 p0 <C 1. The equation of motion for the wavefunction takes the form 



dt 



2m (j) 



V 2 + V B 



0(x, t) + #+(x, £)V(x, t)V(x, t). (35) 



To solve this equation, a mean-field approach is adopted in which 

0(x,t) =^(x,t) + 5V(x,t), 



(36) 



where ^ = (-0) is the phion condensate wave function and Sip is the thermal quantum 
fluctuation around the mean value. The quantity ip is an "order parameter" that 
describes a broken gauge symmetry; the Abelian gauge symmetry occurs in the 
classical MOG action for a massless phion field, m^ = j^j- Moreover, we have that 
(Sip) = and for small Sip the temperature is assumed to be T ~ 0. 

The basic equation for the condensate phion system is the Gross-Pitaevskii equa- 
tion ma EH EH: 



ih 



^(x,t) 

dt 



2m (f> 



V 2 + K grav (x) + A^#(x,t)p 



iP{x,t). 



A stationary solution for the phion wave function can be obtained from 

ip(x,t) = exp(-ifi chcm t/h)ip(x). 
This leads to the time-independent equation 
ti 2 



2m, 



-V z + Vo 



gravV, 



+ N^iPW 



ip(x) = fi chem ip(x). 



(37) 



(3* 



(39) 



The phion condensate gas before recombination can become unstable for a neg- 
ative scattering length, a < 0, corresponding to an attractive inter-phion potential. 
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A stable local minimum for the energy exists only up to a maximum number of 
phion condensates, and above this number the kinetic energy can no longer stabilize 
the condensate gas against collapse. This instability describes the early formation 
of proto-galaxies before recombination and growth of structure after recombination. 
In practise, the energy per phion condensate has three contributions: the kinetic 
energy, the phion potential energy and the potential energy due to gravity. An 
application of the virial theorem to the proto-galaxy condensate gas gives 

3 

£kin = 2^S rav ' (40) 

where S^ n , and £ gra v denote the kinetic energy, the phion potential energy and the 
gravitational potential energy, respectively. As the number of phion condensates 
increases, a critical point is reaches when the condensate gas collapses and clumping 
can begin and form a proto-galaxy. The repulsive phion-phion forces can take over 
at a certain stage of the gravitational collapse of the phion condensate gas for a 
scattering length a > 0, preventing the formation of numerous unobserved compact 
sub-structures in N-body simulations. 

The critical chemical potential for a transition is given by 

8irah 2 . ,„ . s 

/icchcm = ^3"C(3/2), (41) 

where a is the scattering length. For a galaxy-sized object, we can choose a ~ Ac ~ 
14kpc and ~ 10~ 24 eV/c 2 to give 

^^-lO- 14 *- 1 , (42) 

corresponding to a period t p ~ 10 14 s which is roughly the time it takes for a light 
signal to cross a galaxy, Rg/c ~ 10 12 s. 



5 Phion BEC and Evolution of the Universe 

In the late-time universe, the non-relativistic phion matter locally in bound systems 
such as stars, galaxies and clusters of galaxies is dominated by baryon matter and 
neutral hydrogen and helium gases. The quantum non-relativistic, phion condensate 
is subjected to only the weak fifth force with baryon matter, so that as the universe 
expands there is little or no decoherence of the phion condensate gas. However, 
during the evolution of galaxies ordinary baryonic matter begins to dominate inside 
the cores of galaxies and traces light. This predicts that cold dark matter halos 
do not exist around the cores of galaxies or in clusters of galaxies. On the other 
hand, at inter-galactic cosmological scales the neutral BEC phions dominate the 
density of matter with Qb < Q^. The spontaneous symmetry breaking leading to 
the phion BEC is relaxed inside the cores of galaxies, so that the effective mass 
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of the phions is significantly reduced, ~ 10" 24 eV, corresponding to a Compton 
wavelength X c = — h/m^c ~ 14 kpc. 

Locally, inside the bound galaxy systems well after decoupling and the forma- 
tion of galaxies, the spherically symmetric, static solution of the classical MOG field 
equations, together with the variation and renormalization of Newton's gravitational 
constant G, leads at late times to the effective modification of Newtonian accelera- 
tion for weak fields in galaxies [2 . The expression for the modified acceleration is 
obtained from the static spherically symmetric solution of the equation valid outside 
the matter source 

0o + \<t>'o ~ = 0, (43) 

where ' denotes differentiation with respect to the radial coordinate r and 0o denotes 
the time component of the phion field M . The solution to (|4*H|) is the Yukawa 
potential and the modified acceleration law takes the form: 

where 

K = G n ^M^Mq. (45) 

Here, Gn denotes Newton's gravitational constant, M is the total mass of the galaxy 
and Mq is a constant that measures the strength of the coupling of the skew field 
B^ u to matter, expressed in units of the square-root of a mass. We choose 

G OD = G N (l + a) J (46) 

where 

-if- w 

We obtain the expression for the effective modified Newtonian acceleration law: 

a ( r ) = + -exp(-^r)(l + ^r))]. (48) 

We can express the acceleration in terms of an effective Newtonian law in the form: 

a(r) = - G "<;fM, (49) 

where 

Geff(r) = G N [1 + a(l - exp(-/j^r)(l + ^r))}. (50) 

We note that the integration constant K, a and G e g depend on the inverse square 
root of the total mass M. This non-linear dependence on the mass M cannot be 
interpreted in terms of a classical collection of point masses. However, our classical 
modification of Newton's acceleration law is a consequence of an effective classical 
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limit of the non-linear Schrodinger equation for the phion BEC. Let us consider an 
application of Ehrenfest's theorem: 

^(x) = (p/m), ^(p/m) = (a). (51) 

The Hamiltonian operator obtained from Eq.(|37J) is given by 

h 2 

H = -t^V 2 + K grav (x) + W(x), (52) 
where for a static potential 



W(x) = A^(x)| 2 , (53) 



and — ihV = p. We get 



j f (p) = j:([H, p]> = -((VF grav (x)) + (W(x)», (54) 

where — (VV^ rav ) is the Newtonian force law and — (W) is the quantum, BEC phion 
modification of the Newtonian force law corresponding to the modified acceleration 
law in (gU). 

The effective modified acceleration law (}4*%j) fits remarkably well the rotation 
curve data of galaxies and the mass profiles of X-ray clusters [3J E] . This avoids the 
problem of standard cold dark matter models that predict density profiles within 
galactic cores that are too sharp and "cuspy" jlU]- The LSB galaxies have a small 
contribution from the baryonic mass component, so they are efficient tracers of 
rotation curves, and can help to distinguish dark matter halo models from the mod- 
ified Newtonian acceleration law in MOG. The velocity rotation curves predicted by 
MOG fit well the inner cores of dwarf galaxies and LSBs in term of one parameter, 
the mass-to-light ratio < M/L >, using photometric data 

The MOG acceleration formula predicts that for satellite galaxies the rotational 
velocities and density profiles of the satellites have a Newtonian-Kepler behavior. 



6 Weak Field Equations in Cosmology 

In a homogeneous universe with small perturbations, the metric is given in a spatially 
flat universe in the conformal Newtonian-Poisson gauge by 

ds 2 = a 2 [(l + 2<$>)drf - (1 - 2<$>)5 ik dx i dx% (55) 

where i] denotes the conformal time and $Cl. 

For weak gravitational fields the Poisson equation for the gravitational potential 
$ takes the form: 

V 2 $(x, t) + a 2 (t)Y(x, t) = 47ra 2 (t)G(x, t)5p(x, t), (56) 
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where Y denotes the contributions obtained from (JHJ), the factor a 2 (t) is inserted 
to account for the difference between the FRW co-moving coordinate vector x and 
the coordinate vector a(t)x that measures the proper distances at time t, and Sp 
denotes the size of the density fluctuation. If we assume that for large cosmological 
scales, we can ignore the Y contributions and that G = G rcn = constant, then we 
obtain p2|: 

f 1 

$(x, t) = -4vrG rcn t 2/3 a 2 (t) / d 3 q— exp(iq • x)iV q , (57) 

J q 

where iV q = 5 q /t 2 / 3 , <5 q = 5p q /p q . The co-moving wave number vector is denoted 
by q, which is related to the physical wave number at last scattering k = q/a(t LS ) 
where t^s denotes the time of last scattering. 

The equation for the vector field M for weak fields is given by [2j: 



(M) ~ a 2 (t)/4^(x,t) - a\t) d ^^ = a\t)-J^,t), (OS 



where we have assumed that the coupling constant to and the mass /i^ have their 
renormalized values, and J M = (q, J 1 ) (i = 1, 2, 3) is a matter current. 



7 Angular Distance and Cosmological Parameters 

Before recombination the photons couple strongly to baryons, and after recombina- 
tion hydrogen gas becomes neutral and the photons no longer interact with baryons. 
The free photons move along null geodesies and are described by the distribution 
function /: 

dN = f(x\p h T])d 3 xd 3 p, (59) 

where dN denotes the number of photons at time 7]. In the absence of scattering, 
the distribution function / obeys the collisionless Boltzmann equation 

df dx i df dp i df 

"TT - + ~ i o — 7 + ~i n 7 — 0- (60) 

or] at] ox 1 drj op 1 
The geodesic equations for the photons are 

-d7 = P > -d^=2^ PP - (61) 
The Boltzmann equation takes the form 

df df d§ df 

where rii = —pi/p is a vector that determines the direction on the sky from which 
the photons reach an observer. 
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Under general assumptions, the fractional variation from the mean of the CMB 
temperature observed in the direction n takes the form [TH] : 



AT. 
~T~ 



n 



d 3 k5pk exp^d^n ■ k) 



T(k) + in ■ kg(k) 



(63) 



where d,A is the angular diameter distance from the surface of last scattering, and 
fc 2 5pk is proportional to the Fourier transform of the fractional total energy pertur- 
bation. The form factor J-'(k) arises from the Sachs- Wolfe effect and the intrinsic 
temperature fluctuations, while Q{k) arises from the Doppler effect. 

The angular diameter distance of the surface of last scattering is given by 



(L 



1 



ftc /2 # (l 



sinh 



dx 



l/(l+*i) (f^AX 4 + Q k X 2 + Vtrr.x) 1 / 2 



where zl ~ 1100 is the redshift at the surface of last scattering, — 1 — 
is the curvature parameter, and fl m and denote 



87lG vcn p r , 

m 2 



A 
3H 2 ' 



(64) 
(65) 



Here, G ren is the renormalized value of the gravitational constant determined at 
some time t before the time of recombination t ~ t r : 



Cren — Gq{\ + Z), 



(66) 



where Go ~ Gn is the "bare" Newtonian value of the gravitational constant. In our 
modified gravitational theory, the density of matter is given by 



Qb + ^6 + ^£ 



(67) 



where fl^ and Q$ are the densities of matter associated with the phion BEC density 
and the scalar fields, G, u and p, respectively. 



8 Baryon-Radiation Fluid and Phion BEC Fluid 
Before Recombination 

The fluctuations of the background radiation depend upon the radiation energy 
density fluctuations, 5p 7 /p 7 , and the gravitational potential $. The medium con- 
sists of two components, the coupled baryon-radiation fluid and the non-relativistic 
phion BEC condensate fluid and scalar fields fluid. The light phion condensates 
are electrically neutral and do not couple to photons and couple to baryons with a 
gravitational strength. We shall assume that the phion BEC condensate density 
fluctuations 5$ = Sp^/p^ dominate i.e., 5^ 3> 5b and 8$ 3> 5s, where 5b = 5pb/pb and 
5s = 5ps/ps denote the baryon and scalar field fluctuations of the neutral scalar 
fields G,u and /i, respectively. 
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For an imperfect fluid with energy density p and pressure p, the conservation law 
for the standard energy- momentum tensor is VqT ^ = 0, which leads to the first 
order perturbation equation for the components V a T a = [To] : 

Sp' + 3H(5p + 5p) - 3(p + p)$' + a(p + p)d i u i = 0, (68) 

where 'denotes differentiation with respect to rj, u l denotes the spatial velocity for 
i = 1,2,3. Another equation is obtained from the components of conservation, 
V a T a { = 0: 

- 4 {a\p + p)d i u i ]' - ^isAdiU 1 + A5p + (p + p) A$ = 0, (69) 

where ?7 v i S is the shear viscosity. These two equations hold separately for the baryon- 
radiation fluid and the neutral phion BEC fluid. 

For the phion condensate density p^, the classical pressure p and the shear vis- 
cosity rj vis are equal to zero. From (|S8j) and p^a 3 = const., we obtain 

(fy - 3$)' + adiU 1 = 0. (70) 

The dissipation caused by photon pressure and shear viscosity ?7 v i S cannot be ne- 
glected for the baryon-radiation fluid and leads to Silk damping erasure of the 
baryon-radiation perturbations at recombination. 

For non-relativistic baryons, the energy conservation law, V a T a o = 0, reduces 
to a conservation of total baryon number and we obtain from 



(6 b - 3$)' + adm* = 0. (71) 

The equation for the perturbations in the radiation component, <5 7 = <5p 7 /p 7 , is 
given by 

(5 1 - 4$)' + -adiu* = 0. (72) 
3 

The tightly coupled baryon-photon fluid moves with a single velocity, and we get 



Ss 3 
4 



— = ~5 1 — d~b = const., (73) 



where 5s /s denotes the fractional entropy fluctuations in the fluid. Adiabatic per- 
turbations give 5s = and 5b = (3/4)<5 7 . 

The speed of sound in the baryon-photon fluid is given by 

5p 5p 1 + 5pb 3 ( ~^ 4 p 7 ) ^ ^ 

The shear viscosity r] vis is given by 

4 

r/vis = y^P^/, (75) 
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where r 7 is the photon mean-free time. 

The dissipation scale for co-moving wave-number is 

We see from this equation that the viscosity damps the perturbations for co-moving 
scales A < l/kr>. For c 2 = 1/3 and assuming instantaneous recombination, the 
dissipation scale is given by [TH] : 

~r~ ~ 0.6(^) 1 / 4 (fi 6 )- 1 / 2 (^ r )- 3 / 4 , (77) 
k D r) r 

where rj r and z r denote the time and redshift at recombination, respectively. 



9 Correlation Function and Multipoles 

The correlation function for the temperature differences over the sky for a given 
angle 9 is 

= /^(nO^iiA (78) 

where < .. > denotes the average over all ni and 112 with ni • ri2 = cos# and the 
monopole and dipole contributions have been subtracted. We write the correlation 
function in the form 

1 00 

C ( 6 ) = ^E( 2/ + 1 )™(cos^), (79) 



^ 1=2 



where 



S k(Vr)\ ■ n \ 35 k(Vr) dji(kri ) 



0) 



a^/^(^) + ^>(^ Ak d{kvo) 

and t]q denotes the present conformal time. The P;(cos#) and ji{krj) are the Legen- 
dre polynomials and spherical Bessel functions, respectively. In terms of spherical 
harmonics G\ = (|a; m | 2 ) and for 9 < 1 the dominant contribution is for I ~ 1/9. 

In the standard calculations of the CMB power spectrum, the results depend on 
the various cosmological parameters. The generic inflation prediction is, |$|A; 3 | = 
Bk na ~ x , with 1 — n s ~ 0.03 — 0.08. The amplitude B is not predicted and has to be 
fitted to the observations and we assume a flat spectrum n s — 1. The total density 
of matter is given in our modified gravity theory by 

n m = n b + n u + n^ + n s , (8i) 

where 

Oi = ^I^Ei. (82) 

Here, % denotes the contributions from the baryons, neutrinos, phions and scalar 
field components. We have assumed that the gravitational constant has reached a 
constant renormalized value (1661) before recombination. 
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10 Calculation of the Acoustic Oscillation Spec- 
trum 

Mukhanov 1.6] has obtained an analytical solution to the amplitude of fluctuations 
for I > 1: 

l(l + l)Ci (O + N). (83) 

7T 

Here, O denotes the oscillating part of the spectrum, while the non-oscillating con- 
tribution can be written as the sum of three parts 

N = Ni + N 2 + N 3 . (84) 

The oscillating contributions can be calculated from the formula 



O 




A i cos( lr v + ^ j + A 2 cos\2lr p + ^ 



exp(-(/// s ) 2 ), (85) 



where and r p are parameters that determine predominantly the heights and po- 
sitions of the peaks, respectively. The A\ and A 2 are constant coefficients given in 
the range 100 <l < 1200 for Vl b < fi by 

A^0. 1 C ( ' P ^> 2 1/ : 4 - 3 > exp(l( C --/ 7 y). (86) 

•5 + 0.36/ 

;i + 1/2 



A 2 „ 0, 4 (£™P!, (87) 



where 

P = ln( " 



200(^) 



1/2 



and / is given by the ratio 



Vx I /^a\ 1/6 / f y dx V 1 1 . . 

3 7T /_ T-TZXTe) (89) 



Vo zl' 2 \%J Wo (sinhx) 2 ^; ■ 

Here, rj x and z x denote a time and a redshift in the range f] > rj x > rj r when 
radiation can be neglected and y = sinn -1 ^^/^) 1 ^ 2 - To determine r] x /r] , we use 
the exact solution for a flat dust-dominated universe with a cosmological constant 
A: 

a(t) = a (smh^-)H t) , (90) 

where ao and H denote the present values of a and the Hubble parameter H. A 
numerical fitting formula gives [TH] : 

P ~H 2 00(^W/ 2 )' r P=^J d ^(v). (91) 
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Moreover, 



where 



{s ss- 1 = ^ 



3 fpb_ 

4 vp 7 



c s {v) 



1 

71 



1 + e 



(v) \r 1/2 



a(r] r 



For the matter-radiation universe: 

a{rj) = a 

where for radiation-matter equality z = z, 

z, 



+ 2 



+ 2 



(92) 
(93) 

(94) 
(95) 



and r] eq = /^(v^ — 1) follows from a = a(rj eq ). 

The speed of sound at recombination depends only on the baryon density, de- 
termining the deviation from its value in a purely ultra-relativistic medium, and it 
can be expressed as c 2 = 1/3(1 + £). 

The If and l s in (f8T)|) denote the finite thickness and Silk damping scales, respec- 
tively, given by 



I 2 



1 '»V. a 



2a 2 V77 



where 



2(a 2 + l/{k DVr f) 

1/2 



VoV 



a ~ 1.49 x 10 
For the non-oscillating parts, we have 



H (1 + ^ 

z r 



JP -0.22(1/1 A - 3 -2.6) 2 , „„ N9x 

^1 ~ °- 063 C- , , r}Jrnn m, exp(-(/// / ) 2 ), 



H-0.65CZ/Z/) 1 - 4 

0.037 P-0.22(/// s ) a3 + 1.7) 2 
(1 +0 1/2 1 + 0.65(/// / ) 1 - 4 
0.033 P-0.5(/// s ) a55 + 2.2) 2 



exp (-(/// 



, 2 ) 



exp(-(//Z s 



(i + 3/2 i + 2(//g 2 

Mukhanov's formula for the oscillating spectrum is given by 



C{t) = 



1(1 + l)d 100 
[l(l + l)Q} lowl = T 



(O + N), 



(96) 
(97) 

(98) 

(99) 
(100) 

(101) 



where we have normalized the power spectrum by using for a flat spectrum with a 
constant amplitude B: 

9B 

w+mui = Yo^- (102) 
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We adopt the parameters 



Q b ~ 0.04, fi ~ 0.3, Q A ~ 0.7. (103) 

We shall not attempt to separate the degeneracy of the parameter h in the Hubble 
expansion parameter: H = 100 /ikm/Mpc/sec from the parameters f^, ^ and Q A , 
for this requires very accurate CMB power spectrum data, which is not available 
at this time. We adopt in the following the value h ~ 0.71. Moreover, we shall 
not attempt accurate estimates of the parameters and r p . It is important to 
emphasize that the positions and heights of the peaks depend sensitively on the 
values of the cosmological parameters fib, Q<f> arid Indeed, a fit to the CMB 
WMAP, Archeops and Boomerang data requires that fit, < fl m ~ fl^ < Q\ and that 
the universe is spatially flat with Q b + fl + fl A = 1 [TT1 IT21 IT3ll4"Tl 141] . 
The oscillating power spectrum is obtained from the formula: 

C(l) = ^[119.11(0.053((ln(0.0096/) -0.78) 2 -4.3) exp(2.13 x 10~ 7 / 2 ) cos(0.01/ + j) 

+0.14(0.5 + 0.36 ln(0.0096Z)) 2 cos(0.021/ + j)) exp(-8.26 x 10~ 7 / 2 ] 

0.25(ln(0.0096/) - 0.024/ - 3 - 2.6) 2 exp(-4.01 x 10~ 7 l 2 ) 
+ 1 + 5.59 x 10" 5 / 1 - 4 

0.28(ln(0.096Z) - 0.027Z 03 + 1.7) 2 exp(-8.26 x 10~ 7 / 2 ) 
H 1 + 5.59 x 10" 5 / 1 - 4 

0.18((ln(0.0096Z) - 0.011/ - 55 + 2.2) 2 exp(-8.26 x 10~ 7 Z 2 ) 
+ 1 + 3.55 x 10- 6 / 2 ' ( } 

The fluctuation spectrum determined by the analytical formula, Eq. (jl04|) . is 
displayed in Fig. 1 for the choice of cosmological parameters given in (|103|) . 

Power spectrum 



6- 




I 

200 400 600 800 1000 1200 

1 
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Figure 1. The solid line shows the result of the calculation of the power spectrum 
acoustical oscillations: C(l), and the os correspond to the WMAP, Archeops and 
Boomerang data in units fiK 2 x 10~ 3 [TT| IT^l IT3]. 

Our predictions for the CMB power spectrum for large angular scales correspond- 
ing to / < 100 will involve the integrated Sachs- Wolfe contributions obtained from 
the gravitational potential. These predictions should be calculated using a computer 
code and will be investigated in a future article. 

11 Conclusions 

We have demonstrated that a modified gravity theory j2] can lead to a satisfactory fit 
to the acoustical oscillation spectrum obtained in the WMAP data ^Tj by employing 
the analytical formula for the fluctuation spectrum derived by Mukhanov [TSj. The 
vector field M (phion) and the scalar fields G, to and fi occur as new degrees of 
freedom in the action of the MOG theory. The imperfect fluid before recombination 
consists of two components, the baryon-photon fluid and the neutral, non-relativistic 
phion BEC fluid. The latter low-temperature, quantum state condensate fluid has 
zero classical pressure and shear viscosity. 

The analytical formula of Mukhanov reproduces approximately the correct 
peak locations and peak heights in the power spectrum data for the parameters: 
tt b ~ 0.04, tt m ~ ~ 0.3, tt A ~ 0.7 and H = 71 km/Mpc/sec. The varying 
constants G, uj and fi = 1/ro = 1/X are assumed to have reached their renormalized 
constant values before and shortly after recombination. The analytical power spec- 
trum fluctuations are accurate enough to verify that our MOG gives a reasonable 
description of the CMB spectrum for small angular values. 

In the period before recombination and in cosmology at large inter-galactic scales 
the quantum BEC phion energy density dominates the matter density. The BEC 
fluid has zero classical pressure and zero viscosity, whereby the phion field perturba- 
tions as well as the scalar field perturbations are not subject to the Silk dissipation 
that erases the baryon-photon perturbations. The BEC fluid can describe a uniform 
density distribution that does not couple to radiation and can become unstable due 
to the attractive gravitational field and clump to form the seeds of galaxies well be- 
fore recombination, allowing enough time as the universe expands to form clusters 
and super-clusters of galaxies. 

An important problem to investigate is whether an N-body simulation calcula- 
tion based on our phion BEC scenario can predict the observed large scale galaxy 
surveys. The formation of BEC proto-galaxy structure before and after the epoch 
of recombination and the growth of galaxies and clusters of galaxies at later times 
in the expansion of the universe has to be explained. The BEC condensates are not, 
perhaps, localizable as particles in the sense of heavy WIMPS or axions. 
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For local late-time bound systems such as galaxies and clusters of galaxies the 
symmetry breaking is relaxed and the phion condensates can become ultra-light 
and relativistic. Ordinary baryonic matter and neutral hydrogen and helium gases 
now constitute the dominant form of matter. The extra degree of freedom in MOG 
associated with the vector field <fi^ and the spatial variation of Newton's G modifies 
for late-time bound systems the Newtonian acceleration law for weak gravitational 
fields. A calculation of the classical limit of the non-linear Gross-Pitaevskii equation 
for the condensate wave function leads to an effective, classical modified Newtonian 
acceleration law. The rotational velocity curves of spiral galaxies are flattened, 
because of the altered dynamics of the gravitational field at the outer regions of the 
galaxies and not because of the presence of a dominant dark matter halo. 

The neutral phion BEC fluid dominates the cosmological density at large cos- 
mological scales, and the phion condensate in cosmology has taken on the role of a 
light, cold "dark matter" fluid. This allows for a fitting of the acoustical fluctua- 
tion spectrum at small angular scales. According to our MOG scenario, we do not 
expect in the present universe to detect heavy WIMPS such as the supersymmetric 
neutralinos in accelerators or in underground experiments. The dual role played by 
the phion particle in describing galaxies and the large-scale structure of the universe 
is a generic feature of our MOG theory. 

We have succeeded in fitting in a unified picture a large amount of data over 16 
orders of magnitude in distance scale from Earth to the surface of last scattering 
some 13.7 billion years ago, using our modified gravitational theory. The data fitting 
ranges over four distance scales: the solar system, galaxies, clusters of galaxies and 
the CMB power spectrum data at the surface of last scattering. 
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